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Quasiscaling in the analysis of the yield ratio it /tt+i 
Mathematical structure and estimation of source size 

Takeshi Osada, * Minoru Biyajima ^ and Grzegorz Wilk ^ ^ 
Department of Physics, Faculty of Science, Shinshu University, Matsumoto 390, Japan 

(February 1, 2008) 

Recently we have found that integral of the squared Coulomb wave function describing system 
composed of charged pion and central charged fragment Z(,ff protons, |V'r(r)| , times pion source 
function p{r) (of the size /3), f dr \'ipr{r)\'^ p{r), shows a quasiscaling behavior. This is approximately 
invariant under the following transformation: (/3, Z^g) —* (A/3, XZ^ff); A > 0. We called such behavior 
/3-Z(,ff quasiscaling. We examine this quasiscaling behavior in detail. In particular we provide a semi- 
analytical examination of this behavior and confirm it for the exponential pionic source functions in 
addition to the Gaussian ones and for the production of K mesons as well. When combined with the 
results of the HBT, a result of the yield ratio allows us to estimate the size of the central charged 
fragment (CCF) to be 125 < Z^g < 150 for Pb+Pb collisions at energy 158 GeV/nucleon. From 
our estimation, the baryon number density 0.024 < n^ < 0.036 [1/fm"^] is obtained. 



I. INTRODUCTION 



Recently E866 1 1| and NA44 M collaborations reported a remarkable deviation from unity of the yield ratio n /n^ 



measured at smal 



niT-iTiTr region in high energy heavy-ion collisions. Gyulassy and KaufFmann p] have argued that 

T^ , such deviations are due to the Coulomb interaction between charged n^ and a system of Z^a protons (called a central 

fT^ ' charged fragment - CCF) which are nearly stopped in the central rapidity region in the course of the heavy ion 

C^l [ collision. They have provided theoretical formulae for this effect calculated up to the first order in the ZcSOi- On 

^^ ■ the other hand Shuryak M has also pointed out that the above deviation seen in the high energy nucleus- nucleus 

l^ ' collisions can be regarded as a manifestation of a Coulomb interaction effect and have proposed a theoretical formula 

^^ . using the Gamow factor. See also refs. |^,|| 

^ ; 

Mh, In our previous paper J7| we have derived a new theoretical formula for the observed yield ratio 7r^/7r+ using the 

O , Coulomb wave function. We have found that, when applied to the experimental data, this formula exhibits a valley 

^ , like structure of the minimum x^ values in P-Z^b parameter space. This result suggests existence of the following 

1-^ fi-Zcs quasiscaling in the pion production yield iV" 



,± 



N'^ {rriT-mT,; Xxr]±, Ax/3) w N'' {rriT - m.,:; t]±, /?), (1.1) 



j^ • where ri± — iZcs fJ.a/Pr ^^d N^ is given by 



N''^{mt-m^-,r]±,P) = / d^r p{r) iptiPr,r) 



CXD CXD 
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G'(^±)E E W;^ MV±)AUV±) Hn.m) \ 2p. 1"+". (1.2) 
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Here G{r]±) — 27r77-t/(exp(27r77±) — l) is the Gamow factor and 



^"(^±) = -r7^^^T^' Ip{n,m)^47rldrr"^+"^+'p{r). (1.3) 



r(i?7±) (n^-) 
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The Coulomb wave function ip^{pr,r) of the CCF system of Zcff protons and a charged n^ is a solution of the 
following Schrodinger equation: 



2// r 



^t{v)^Eri^f{r), (1.4) 



and is given explicitly in terms of the confluent hypergeometric function F |9[| : 

^±(p,,r)=^(l + ^r?±)e-''■''±/2e'P-r^(_^r;±,l,i(p,r-p,•r)) . (1.5) 

For the Gaussian source, p{r) — I .J" ) exp I -^^ ) and for the exponential source, p{r) = "g^m cxp {-g-) , we have 
Ip = 2 ^^^)n+m r( "+^+3 ) and If) = i/?"+™r(n + TO + 3), respectively In the next section we demonstrate that 
the quasiscaling behavior is independent of the type of the pionic source used. In Section III we examine in detail 
the quasiscaling behavior of both N"^ and N^ production yields. In particular, we provide both semi-analytical 
and numerical proofs of this behavior. In Section IV we estimate the most probably size Z^s of the CCF and baryon 
number density of it using results of tht HBT analysis. In the last section we present our concluding remarks. 



II. ANALYSIS OF DATA BY NA44 COLLABORATION USING EXPONENTIAL SOURCE FUNCTION 

We analyze data of the yield ratio tt^/tt^ observed in Pb+Pb collisions at 158 GeV/nucleon g] by applying the 
following formula, 

_ ^ N"" {mj,-m^;r]_,(3) 
' 7V-^(to^-to,;77+,/3) ^ ' 

The exponential source function in TV is used in eq.(EJ). In Fig.l we present the result by this source function and 
the previous result in ref. . In both cases we fix values of /3 and estimate Zcff providing the minimal value of x^ 
(by using CERN-MINUIT program). In both cases the corresponding parameters Z^ff and /3 are proportional to each 
other: Z^ff ~ k/3 (where k is a constant depending on the type of source function used). In this way we confirm that 
the quasiscaling behavior is independent of the type of source function used. 



III. EXAMINATION OF THE QUASISCALING BEHAVIOR OF A^'^ AND N^ 
A. Taylor series expansion and finite difference method 

We elucidate the mathematical structure of the quasiscaling behavior seen in Fig. 1. The numerical results presented 



there suggest that the series expansion part of eq.(1.2) with arguments A?7± and A/3 has the same A dependence as 
the 1/G'(A?7±) part, i.e., that 



Six 7y±, A /3)ee^ Y. \ !, A„(A77±)A;,(Ar;±) I(^xp){n,m) \ 2p, 
^ — ' ^ — ' n + TO + 1 

n— m—0 

^constant x l/G{Xri±). (3.1) 



To examine this supposition we expand the 5(A?7+,A/3) (with the Gaussian source function) and 1/G'(A?7+) in the 
Taylor series around Aq {— 1,2) with fixed mr-TO,^(= 50, 100 and 200 MeV) by using the Mathematica software. All 
combinations of indices {n,m) satisfying < n + to, < 90 were considered. The result is presented in Table || [ [lO[ . 
One can observe that coefficients of A expansion in l/G{Xri±) agree reasonably with those emerging from S{Xri^, Xf3) 
(modulus constant multiplying S{Xr]+, A/3).). 

However, it is difficult to evaluate higher order terms of the series in S{Xri^, A/3) using this method [Q. Therefore, 
we also calculate five coefficients Cj {0 < j < 4) of the A-expansion in 5(A?7+,A/3) and 1/G'(A77+) using the finite 
difference method (with its size equal to e = 1/200) defined as follows : 



F{X) ^co + ci SX + C2 SX^ + C3 ^A^ + C4 JA* + 
where ^A = | A — Aq | -^ e and , 
CO = F{Xo), 

l-[F{Xo+e/2)-FiXo-e/2)], 
^ ■[F(Ao+e)-2F(Ao)+F(Ao-£)], 



(3.2) 
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C3 
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C4 = ^ [F(Ao + 2£) - 4 F(Ao + e) + 6 F(Ao) - 4 F(Ao - e) + i^(Ao - 2£)] . 



[F(Ao + 3£/2) - 3 F(Ao + £/2) + 3 F{Xo - £/2) - F{Xo - 3e/2)] , 



Comparing coefficients of A-expansions in S'(A 77, A /3) and l/G(Xr]) (cf. Table ||), we confirm that e c|.(|3.lD approxi- 
mately holds. This in turn means that the A dependences of G{Xri) and S{X 77, A /3) in the l.h.s. of eq.( |l.l| ) effectively 
cancel out and, as the results, the quasiscaling behavior emerges. The same is true for the exponential source function. 



In Tables HI, IV we present results of same examinations of S{X ?7_|_,A /3) and l/G(Ar/+) but for the production 



of K~^ meson instead (Table [II contains results of the the Taylor series expansion and Table IV those for the finite 
difference method). As seen in Tables [II, [V, they confirm that coefficients in l/G{Xrj) agree with those of S{X r], X (3) 
also for the K meson production case. Notice that eq.(3.1) holds well for large m.r-m^(K) for both N^ and N^ 
(this is clearly seen in the first order e of A=2 in Tables I - IV). 



T-TriK) 



We examine dependence of the coefficients of S{Xri+, A/3) and l/G{Xri+) expansions on the size of e (both for the 
Gaussian and for the exponential source functions), cf. Table [v| and |V^(5'(A 77+, A /3) calculated using the Gaussianand 
the exponential source functions are denoted by Sg and S'e, respectively). Expansions are done around Ao(=2) with 
fixed values of 171^,-171.^ or rriT-mK (= 200 MeV in both cases) with e — 1/2, 1/20 and 1/200, respectively. It can 
be observed that coefficients of the A-expansion in l/G'(Ar7+) agree reasonably well with those for S{Xr]+,XP) up to 
£ = 1/2. This means that quasiscaling behavior holds not only on the local scale (i.e., for small £) but also globally 
(up to £ of the order of 1/2). 



B. Numerical examination of the quasiscaling behavior 



To visualize the global structure of the f3-Zcs quasiscaling more clear we examine the A-dependence of the following 
quantity: 



£ 



N^\iX + e/2)rj, iX + e/2)p) - N^\ (X- e/2)rj, (A - £/2)/3) 



(3.3) 



Figure 2 shows the A dependence of the I?^ (A; m^ — ttitt) for both the Gaussian (a) and the exponential (&) source 
functions used for N'^ . Figure 3 shows the same but for the K~^ meson production instead, i.e., the A dependence 
of the D^ (A; ttt-t — ttt-k)- Except for the cases in which either A ^ 1 or mT-mT^f^x) is small, D'^ ^^ ^(A; nij. — 111^^) 
is e ssen tially equal zero. In these figures we can readily observe disappearance of the A dependence of the l.h.s. of 
eq.(lT) and the emergence of the quasiscaling behavior. On the other hand we can also observe that for small values 
of m-r — m^(^x) the quasiscaling breaks down because D{X] m^ — rriT,) ^ there p2|. 



IV. HOW TO ESTIMATE THE SIZE OF THE CCF 

It is tempting to use our results to estimate the size of CCF as given by the effective number of stopped protons, 
Zcff. To this end we need some additional inputs, which can be provided in two ways: Either by adopting empirical 
geometrical picture of the nucleus or by utilizing results obtained from the HBT effect. 



First we use the empirical formula for the radius of a nucleus (of an atomic number A) , 



{rl)^'^^r,A^'\ (4.1) 

where rg = 1.1 ^ 1.2 and, in our case, A « (2.5 ~ 3.0) x Z^s, because of A=207 and Z=82. We use the root mean 
square radius of the pionic source as (rj)^'^ = {47r J (irr'*p(r)}^'^ (which in terms of the source size parameter (3 
used before is given as (r^)^'^ =v3/3 and (rj)^" —2^/2>f3 for the Gaussian and the exponential sources, respectively.). 
Using (rj)^/^ instead of /3, we can rewrite Fig.l. cf. Fig. 4. Dotted line shows there the relation 



^cff ■ 



' {rlf'\ (4.2) 



2.5r3 



which provides us with Z^s = 45 ^ 60 for (r^)^/^ w 6 fm (where both lines cross each other). 

In the second method we first apply the following standard formula with the Gaussian source function to the data 
of Pb+Pb collisions at energy 158GeV/nucleon W%: 

iV±±/7V^G=c(l + Ae-'3'Q'), (4.3) 

(here c is a normalization constant and A is a degree of coherence parameter). As a result we obtain the following 

sets of parameters for tt+tt"*" and tt^tt^: 

(c, A,/5) = (1.00±0.01, 0.55±0.03, 6.02±0.35) with xVNDF=13.0/17 and 

(c, A,/3) = (1.01±0.01, 0.47±0.03, 5.58±0.40) with xVNDF=7.4/17, respectively. Therefore, the pionic source size of 

Pb+Pb collisions is estimated to be equal 12.7 fm< (rj)^/^ < 15.6 fm. Using this interval as the size of the CCF, we 

estimate that 125 < Z^s < 150, which is much larger than the that of CCF obtained from the empirical geometrical 

picture. It means that the CCF represents an expanding source, provided that the estimation by means of HBT effect 

is correct. From our estimations, the baryon number density 

„ , ,.„ / 47r( 15.1 fm)3 ^ ^ ,„. / 47r( 12.5 fm)^ 
2.5 X 150 / — i —<nB < 2.5 x 125 / — ^ — 

0.024<nB < 0.036 [1/fm^] (4.4) 

is obtained. 



V. CONCLUDING REMARKS 

The mathematical structure for the f3-ZcS quasiscaling in eq. (fj) found previously in W is examined in more detail. 
In particular, through the Taylor series expansion (up to third order) we find that coefficients in the A-expansion of 
1/G(A77) are approximately same as those of S{X ri,X (3). Using the finite difference method we can improve our 



examination by extending it to the fourth or fifth orders, with similar result. This means that eq.(3.1) holds indeed 



(albeit only approximately) and that A-dependence of the l.h.s. of eq.(l.l) is cancelled out in the product of G(Ar/) 
and ^(A 77, A /3) leading in consequence to the quasiscafing behavior observed. Finally we estimate the Zes of the CCF 
in two possible ways. From our estimations, the baryon number density 0.024 < ne < 0.036 [ 1/fm ] is obtained. 
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FIG. 1. Valley structure of the minimum x^ values in Zatt-P parameter space for NA44 Collab.[2] data for the Gaussian and 
exponential pionic source functions, respectively. 

FIG. 2. Numerical examination of the Z^n-P quasiscaling in A"^ {mr — m.^; Xxrj, Ax/3) for the Gaussian (a) and exponential 
(6) source functions. 

FIG. 3. The same as in Fig. 2 but for the production of K mesons. 
FIG. 4. Two lines in Fig. 1 (for the Gaussian and the exponential source function) are rewritten by using of the r.m.s. radius 



r: 



2^1/2 rpj^g dotted line is evaluate by eq.(ll) 



TABLE I. Comparison of the coefficients of the A-expansion (the Taylor series) in 1/G'(A?7) and in S{X rj, A /3) 
around the Ao with fixed mT-m,r; 5X = | A — Ao |<C e. The Gaussian source function is used in ^(A r;, A /3). 
In this Table we omit imaginary parts of coefficients in the expansion of ^(A rj, A /3) emerging from the limited com- 
puter ability. For example, in A=l, mT-m,r=50 MeV case, we obtain in reality following result at out computer facility: 
S{\ri, A/3) = 1.6531 + i 1.1956 x 10"^** + (0.8622 + i 1.6171 x lO"^'^) e + (0.3853 + i 1.0199 x 10"^®) e^ + • • ■ . 



A*"^ mx-Jn^ 


Ao = 1 50 MeV 

1/G(A?7) = 1.6531 + 0.8902 S\ + 0.2971 SX"" + ■ 

S{\rj, A/3) = 1.6531 + 0.8622 5\ + 0.3853 S\^ + ■ 




Ao = 1 100 MeV 

1/G(A?7) = 1.4129 + 0.5134 5\ + 0.1180 3\^ + ■ 

S{\r], \f3) = 1.4126 + 0.5274 SX + 0.1217 5X^ + ■ 




Ao = 1 200 MeV 

l/G(Ar?) = 1.2468 + 0.2845 SX + 0.4180 ^A^ + • 

SiXr], A/3) = 1.2467 + 0.2873 5X + 0.3674 SX^ + ■ 




Ao = 2 50 MeV 

l/G(Ar?) = 2.9283 + 1.7687 SX + 0.6241 ^A^ + • 

S{Xri, Xf3) = 2.9272 + 1.7866 5X + 0.6238 SX"" + ■ 





TABLE IL Comparison of the coefficients of the A-expansion (the finite difference method, e — 1/200) in ^(A rj, A /3) around 
the Ao with fixed m-r-mK- The Gaussian source function is used in ^(A r],X /3). 



N" 



mj--mTr 



Ao = 2 


50 MeV 


1/G(A7?) 


= 2.9283 + 1.7687 SX + 0.6242 SX^ + 0.1572 SX^ + 0.0309 (JA"* + ■ ■ ■ 


S{Xn,XP) 


= 2.9283 + 1.7867 <5A + 0.6239 SX'^ + 0.1555 ^A^ + 0.0358 ^A"* + ■ ■ ■ 


Ao = 2 


100 MeV 


1/G(Aj?) 


= 2.0672 + 0.8216 SX + 0.1973 ^A^ + 0.0342 5A^ + 0.0047 SX'^ + ■■■ 


s{Xn,xp) 


= 2.0672 + 0.8245 <5A + 0.1958 SX'^ + 0.0356 SX^ + 0.0037 ^A"* + ■ ■ ■ 


Ao = 2 


200 MoV 


l/G(Ar,) 


= 1.5780 -1- 0.3834 SX + 0.0581 SX^ + 0.0064 5A^ + 0.0006 ^A"* + ■ ■ ■ 


5(A7?,A/3) 


= 1.5780 -1- 0.3837 SX + 0.0579 SX^ + 0.0066 SX^ + 0.0005 ^A"* + ■ ■ ■ 



TABLE IIL The same as in Table I but for the production of K mesons. 



N'- 



rrij^-mK 



Ao = l 

1/G(A77) 

S{Xv,Xf3) 



50 MeV 

= 2.2584 + 2.0291 SX + 1.0916 SX^ + 

= 2.2584 + 2.0459 SX + 1.0777 SX'^ + 



Ao = l 

1/G(A77) 

S{Xv,Xf3) 



100 MeV 

= 1.9074 + 1.3380 SX + 0.5723 ^A^ + 

= 1.9074 -I- 1.2893 SX - 1.9604 SX^ + 



TABLE IV. The same as in Table II but for the production of K mesons. 



A'" mT-rriK 



Ao = 1 100 MeV 

l/G(Ar?) = 1.9074 + 1.3381 S\ + 0.5723 3\^ + 0.1774 5X^ + 0.0432 S\^ + 

g(A?7,A/3) = 1.9074 + 1.3417 5\ + 0.5669 5\'^ + 0.1866 SX^ + 0.0298 ^A'* + 



Ao = 2 50 MeV 

l/G(Ar;) = 5.9708 + 6.2065 <5A + 3.5872 5\'^ + 1.4529 SX^ + 0.4551 SX"^ + 

S{Xr], A/3) = 5.9708 + 6.2079 5X + 3.5864 ^A^ + 1.4531 SX^ + 0.4541 ^A"* + 



Ao = 2 100 MeV 

1/G(A??) = 4.0487 + 3.2420 SX + 1.4789 SX^ + 0.4770 SX^ + 0.1195 SX^ + 

SjXy, A/3) = 4.0487 + 3.2420 SX + 1.4786 ^A^ + 0.4777 SX^ + 0.1189 5X^ + 



Ao = 2 200 MeV 

1/G(A??) = 2.7398 + 1.5454 3X + 0.5123 3X^ + 0.1215 5X^ + 0.0225 SX^ + 

S(A??, A/3) = 2.7398 + 1.5463 SX + 0.5172 5X'^ + 0.1296 SX^ + 0.0022 ^A"* + 



TABLE V. Dependence of A expansion of l/G{Xri), Sg{X ri,X /3) and S'e(A r;, A /3) for tt"*" production (in finite difference 
method) on the size of e (for fixed values of Aq = 2 and tht — Jn^ = 200 MeV.) 

N^ rriT-m^ ; ( /3=2 fm, Z^ff^AS ) 



£=1/2 


200 MeV 


l/G(Ar;) 


= 1.5780 + 0.3838 <5A + 0.0582 ^A^ + 0.0065 SX^ + 0.0006 ^A"* + ■ ■ ■ 


S,{Xv,Xl3) 


= 1.5780 + 0.3841 SX + 0.0580 SX^ + 0.0067 5X^ + 0.0004 ^A"* + ■ ■ ■ 


S.(Ar;,A/3) 


= 1.5780 + 0.3839 SX + 0.0581 ^A^ + 0.0065 5X^ + 0.0005 .JA"* + ■ ■ ■ 


e = 1/20 


200 MeV 


l/G(Ar;) 


= 1.5780 + 0.3834 SX + 0.0581 SX^ + 0.0064 SX^ + 0.0006 ^A"* + ■ ■ ■ 


Sg{Xv,Xf3) 


= 1.5780 + 0.3837 <5A + 0.0579 SX'^ + 0.0066 ^A^ + 0.0005 ^A"* + ■ ■ ■ 


Se(A77,A/3) 


= 1.5780 + 0.3835 <5A + 0.0580 SX^ + 0.0068 5X^ + 0.0039 ^A"* + ■ ■ ■ 


£ = 1/200 


200 MeV 


l/G(Ar;) 


= 1.5780 + 0.3834 SX + 0.0581 SX^ + 0.0064 SX^ + 0.0006 ^A"* + ■ ■ ■ 


S, (At?, A/3) 


= 1.5780 + 0.3837 <5A + 0.0579 SX^ + 0.0066 5A^ + 0.0005 ^A"* + ■ ■ ■ 


Se(Ar?,A/3) 


= 1.5780 + 0.3835 <5A + 0.0579 SX^ + 0.0069 5X^ + 0.0047 ^A"* + ■ ■ ■ 



TABLE VI. The same as in Table V but for the production of K mesons. 

A-^^ mT-rriK ; ( /3=2 fm, Zeff=48 ) 

£ = 1/2 200 MeV 

l/G(Ar?) = 2.7398 + 1.5530 SX + 0.5180 SX'^ + 0.1236 SX^ + 0.0231 ^A"* + 

SgiXrj, XP) = 2.7398 + 1.5541 <5A + 0.5189 SX^ + 0.1219 5X^ + 0.0200 ^A"* + 

5e(A?7, A/3) = 2.7398 + 1.5535 SX + 0.5181 SX^ + 0.1225 SX^ + 0.0222 SX'^ + 



1/20 200 MeV 

1/G(A77) = 2.7398 + 1.5455 SX + 0.5124 SX^ + 0.1215 5A^ + 0.0225 ^A"* + 

Sg{Xrt, A/3) = 2.7398 + 1.5464 SX + 0.5172 SX^ + 0.1292 5X^ + 0.0017 ^A"* + 

Se(Xri, A/3) = 2.7398 + 1.5453 SX + 0.5140 SX^ + 0.1504 5X^ + 0.0160 SX^ + 



£ = 1/200 200 MeV 

l/G(Ar?) = 2.7398 + 1.5454 SX + 0.5123 SX^ + 0.1215 5A^ + 0.0225 ^A"* + 

Sg{Xrj, A/3) = 2.7398 + 1.5463 SX + 0.5172 (JA^ + 0.1296 5X^ + 0.0022 ^A"* + 

Se{Xri, A/3) = 2.7398 + 1.5453 SX + 0.5139 SX^ + 0.1541 5X^ + 0.0130 SX^ + 
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